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ABSTRACT 

We report the first detections of circularly polarized emission at subniillimeter wavelengths from the 
compact radio source and supermassive black hole candidate Sgr A* at a level of 1.2 ±0.3% at 1.3 mm 
wavelength (230 GHz) and 1.6 ± 0.3% at 860 /im (345 GHz) with the same handedness, left circular 
polarization (LCP), as observed at all lower frequencies (1.4-15 GHz). The observations, taken with 
the Subniillimeter Array in multiple epochs, also show simultaneous linear polarization (LP) at both 
wavelengths of about 6%. These properties differ sharply from those at wavelengths longer than 1 cm 
(frequencies below 30 GHz), where weak circular polarization (CP) (~ 0.5%) dominates over LP, 
which is not detected at similar fractional limits. We describe an extensive set of tests to ensure the 
accuracy of our measurements. We find no circular polarization (CP) in any other source, including 
the bright quasar 1924-292, which traces the same path on the sky as Sgr A* and therefore should 
be subject to identical systematic errors originating in the instrument frame. Since a relativistic 
synchrotron plasma is expected to produce little CP, the observed CP is probably generated close to 
the event horizon by the Faraday conversion process. We use a simple approximation to show that the 
phase shift associated with Faraday conversion can be nearly independent of frequency, a sufficient 
condition to make the handedness of CP independent of frequency. Because the size of the t — 1 
surface changes by more than an order of magnitude between 1.4 and 345 GHz, the magnetic field 
must be coherent over such scales to consistently produce LCP. To improve our understanding of the 
environment of SgrA* critical future measurements include determining whether the Faraday rotation 
deviates from a dependence in wavelength and whether the circular and linear components of the 
flux density are correlated. 

Subject headings: black hole physics — Galaxy: center — plasmas — polarization — subniillimeter 
— techniques: interferometric 



1. INTRODUCTION 

The Galactic Center source Sagittarius A* (Sgr A*), 
the nearest supermassive black hole (SMBH), is ex- 
tremely underluminous for its mass ('^ 4 x 10^ M©, 
Chez et al.|2008 Gillessen et al.|2009 ), radiating at only 
10~^LEdd- Theoretical models have focused their efforts 
on explaining this faintness by invoking diverse physical 
mechanisms that result in r adiatively inefficient accretion 
and/or ou tflow process es ([Falcke e t al. 1993; Narayan 
FYilfT994; BlandfoTd^fe Begelman 1999; Quataert ^ 
Falcke & Markoff,2000^ Yuan et al.. 20021 



et al.||2008[ [Fish et al.||2011 ). The size- wavelength rela- 
tion, and even the observed interferometric visibility of 
the emissio n, can b e matched by several model s ( Markoff] 
et al |2007[ 'Moscibr odzka et al.|2009l [Huang et al. 2009Id 



BrodericTet aL|201l|). Multiwavelengt h vari ability dMar 
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Yusef-Zadeh et ah 
combination with 



Jruzinov^OOO ; 



2008[ IKunneriath 
m comDmation witn other observables, 
which may also reveal information about sourc e struc- 
ture, can a, l so be replicated in multiple schemes ( Dexter 
et al.||2009| [Maitra et aT1|2009[ ) 



I'hese models a dequately reproduce~the quiescenl; spec 
Sgr A* dNarayan et a l.'1998i [Markoff et alT200T 
; al.||2001t lYuan et ai. 2003), although the spec 



trum of 

Melia et al.||2001[ [Yuan et~a l. 2003), altliougli tlie spec 
trum alone does not discriminate between them. 

In subsequent years, new observations have provided 
new types of constraints on the models. Very long base- 
line interferometry (VLBI) has measured the wavelength- 
dependent size of Sgr A*, detecting structure on event- 
horizon s cales in the highest frequency /resolution exper- 
iments ( )Shen et al. 2005 Bower et al.||2006| iDoeleman 



I'he polarization spectrum and its variability present 
additional, rich information about source structure. 
Strong theoretical constraints can be derived by consid- 
ering the effects of polarized radiative transfer (PRT) in 
the models. Perhaps the simplest and most common con- 
sequence of PRT through a plasma is Faraday rotation of 
linear polarization (LP). In recent years, multifrequency 
millimeter and submillimeter polarimetry ha s allowed 
measurements of Faraday rotation in Sgr A*(Marrone 
et al |2006a[|2007[ ). These data have pr 
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the accretion rate in the inner regions of the flow (|Mac- 

eraI1|2007| ). 
lasma can b( 
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A synchrotron-emitting plasma can be expected to pro- 
duce signiflcant LP fractions at frequencies near the spec- 
tral peak. At the same time, these relativistic plas- 
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TABLE 1 

Radio-to-Submillimeter Linear and Circular Polarization of Sgr A* 



Frequency (GHz) 


Fractional LP (%)=» 


Fractional Stokes V (%)'' 


Reference 


1.4 
4.8 
8.4 

15 




—0.21 ± 0.10 (13) 
-0.33 ±0.07 (13) 
—0.32 ± 0.08 (13) 
—0.62 ± 0.26 (12) 


|Bower et al.|(|2002|| 

„ 


4.8 




— 37 + 04 


Sault & Macquart (19991 


4.8 
8.4 


< 0.08 

< 0.17 




jBower et al.| Ijl999a|) 


4.8 

8.4 




— 36 + 05 (3) 
—0.26 ± 0.06 (3) 


|Bower et ai.| l|1999b|| 


22 
43 
86 


< 0.2 

< 0.4 

< 1.0 




|Bower et al.|l|l999c|| 


82.8 
86.3 


2.1 ± 0.4 
0.8 ±0.5 




|Macquart et al.|||2006|l 


100 




< 2.0 


'I'suboi ct al. (2003) 


112 


< 3.6 


< 3.6 


Bower ct aL (2001J 


216 
230 


9.1 ± 2.2 
10.0 ± 2.5 




|Bower et al.|i2UUb| 


230 


7.2 ±0.6 


^ 2 


•Bower ct al. (2003) 


230 


5.9 ± 1.6 




JViarronc ct al. (2U07I 


340 


6.4 ± 2.0 




Marronc ct al. (^UU6a) 


150 
225 
350 
400 


22^9^ 




Aitkcn ct al. ( 2000 ) 
■ . . 1 

)) 



^ Uncertainties can correspond to both systematic err ors and time- va r iability dispe rsion. For example, 
data corresponding to the multi-epoch observations of |Marrone et al.| j2006a| |2007[ l are presented with 
errors corresponding to the standard deviation of the sample. 

^ The errors shown are the standard deviation (rms) of the sampl e, not the standard deviation of the 
mean. The number of samples is listed in parentheses. For the [Bower et al. data, we show 

the results for the 1999 VLA measurements only because of the advanced calibration technique used. 
Howev er, they are consistent with the archival VLA and the ATCA data also reported by [Bower et ah] 
j2002( . Upper hmits are shown at the 2-cr leveL 



mas are expected to produce very small amounts of cir- 
cu lar polarization (CP) by intrinsic emission of Stokes 
Ff La ndau k Lifshi tz'WzK^ or through various radiative 
transfer ettects (e.g., Jones & Odelliigrn pelrose,.199'n 
Ruszkowski fc Begelman||2UUlj| |Shcherbakov||lj008| ~~5j 
the polarization arises very near to the black hole, a full 
picture of the emergent polarization state can diagnose 
both the inner regions and the intervening propagation 
medium. 

A compilation of published polarization measurements 
is given in Table [l] and shown in Figure [l] No LP has 
been observed at 8 GHz or below at limits of 0.2% or 
less dBower et a r| | 1999a[), The LP of Sgr A* was first 
detected by |Aitken et af] ( |2000| ) above 100 GHz. These 
measurements were obtained with large beams (7"-22") 
and required background subtraction, leaving some un- 
certainty about residual polarization contamination from 
the surrounding dust emission. Subsequent interferom- 
eter observations at millimeter and submillimeter wave- 
lengths have shown polarizatio n at the level of 1- 10%, 



which va ries in position angle ( Bower et al. 2005 1 and 
fraction ( Marrone et al. 2006a ) with timescales compa- 
rable t o those ot previously ob served total intensity vari- 
ations (Marrone et al. 2006b). The variability may be 
intrinsic to the source or due to propagation effects, but 
the short timescales involved suggest that processes very 
close to the SMBH are responsible. The LP was mea- 
sured si multaneously at multipl e frequencies for the first 
time by Marrone et al. (2007). The inferred rotation 



measure (RM) indicates that Sgr A* is extremely under- 
fed, with an accretion rate of ^ 10~* M© yr~^. 



CP f rom Sgr A* was first detecte d by Bower et al. 
(|1999bj). Later, iBower et al.l ([20028 reported spectral 



measurements between 1.4 and 15 GHz as well as time 
variability of CP (see lower panel in Figurejl]). They con- 
cluded that the time-averaged CP spectrum is approxi- 
mately i/0-5±0-2^ with persistent variability that increases 
with frequency. They also noted that the sense of CP 
for all available measurements (about 100 measurements 
from 1981 to 1999, made with both the VLA and ATCA 
at 1.4, 4.8, and 15 GHz) was exclusively LCP, indica- 
tive of a long-term stability in magnetic field configura- 
tion. At these frequencies, LP is not detected down to 
instrumental limits of 0.1% while CP is persistently de- 
tected at levels of a fraction of a percent with a degree of 
variability of the same order of magnitude. This is sub- 
stantially different from what is seen in high-luminosity 
active galactic nuclei ( AGN) , where LP always dominates 
CP. Above roughly 100 GHz, LP becomes dominant in 
Sgr A*, with CP undetectable at ~1% sensitivi ty. In 



four epochs of 230 GHz observations, Bower et al. (20031 
found a 2-cr CP signal (3%) on one day, with an average 
measurement of 1%± 1%. No CP was observed at a level 



-0.5 ±0 3^ 



of 0.5% in the 100 GHz data reported by |Tsuboi et al. 
( |2003| . [Marro ne et"aL] ([2006a[) measured V at 345 GHz 
and obtained 



However, because of system- 
atic calibration uncertainties, they reported this result 
as an upper limit of 1.5%. 



CIRCULAR POLARIZATION OF SCR A* 



3 



c 
o 

o 

[In 



100 



^ 10 



Bower et al (1999a,c; 2001) 
Macquart et al (2006) 
Aitken et al (2000) 
Bower et al (2005) 
Bower et al (2003) 
Marrone et al (2007) 
Marrone et al (2006a) 



0.1 



X 
-B- 



Oh 

u 

a 
_o 

o 



10 



0.1 



Bower etal (2001) 
Tsuboi et al (2003) 
Bower et al (2003) 
Marrone et al (2006a) 
Sault& Macquart (1999) 
Bower etal (1999b) 
Bower et al (2002) 



10 



100 



Frequency [GHz] 



Fig. 1. — Published measurements of the fractional linear and circular polarization toward Sgr A* as a function of frequency. 



CP is also observed in a variety of radio sources, includ- 
ing pulsars and AGN. Examples are 3C273 and 3C279 
(Homan & Wardlc 1999). Some models seeking to ex- 
plain the millimeter and submillimeter LP have also pre- 
dicted CP at these high frequencies du e to the conversion 



of LP to CP in a turbulent jet (e.g., Becker t fc Falcke 
20021. In these models, in addition to the stochasticity 
of the magnetic field — which appears to play a crucial 
role in building up CP by propagation effects — the heli- 
cal geometry of jets might be important in high levels of 
LP, above 100 GHz. 

Coupled mechanisms to produce both LP and CP in 
relativistic outflows ha ve been st udied in detail by many 
autho rs (e.g., Ruszkowski & Begclman 2002; Bec kert fc 
Falcke 2002; Beckert 2003; Huang et al. 2008; Homan 
et al.|i2009: Shcherbakov et 



al. 



. 2010^ . If CP is pro- 
duced pi^ominannjr'Bjr^ropa^ effects, multifre- 
quency measurements of RM and CP can provide im- 
portant clues about the magnetic field structure of the 
plasma surrounding the SMBH. Similarly, simultaneous 
measurements of LP and CP variability can determine 
whether the intraday variability in Sgr A* is due to in- 
trinsic variations in the central engine or to variations 
in the outer layers of the accretion flow. In this way, it 
should be possible to conclusively infer the prese nce or 



absence of a "Faraday screen" in front of Sgr A* (Mac- 



quart et al.||2006| [Marrone et al.| 2007| ). 



in Section [2] we describe several epochs of polarimetric 
observations of Sgr A*, Section [s] reports the observed 
polarization, and Section [4] explores the tests of system- 
atic errors in the observations that might give rise to 
false CP. Section [5] reviews PRT and the astrophysical 



TABLE 2 
Observation Epochs 



Date 


Main 


Freq.^ 


SMA 


''"225 GHz 




Target 


[GHz] 


Config.*^ 




2005 June 6 


Sgr A* 


343.0 


CN 


0.055 


2006 July 17 


Sgr A* 


226.9 


VEX 


0.05-0.08 


2007 March 31 


Sgr A* 


226.9 


CN 


0.055 


2008 May 30 


1924-292 


226.9 


CN 


0.08 



^ Frequency of the local oscillator. Upper and lower sidebands 
are centered 5 GHz above and below this frequency, respec- 
tively. 

^ Array configurations irrclude "Compact North" (CN) and 
"Very Extended" (VEX). 

mechanisms for generating CP, and their applicability to 
Sgr A*. 

2. OBSERVATIONS 

The data presented in this work consist of three po- 
larimetry tracks (Table |2| taken with the Submillimeter 
Array (SMA) H The gener a l char acter istics of th e SMA 
are described by |Blundell| ( |2004| a nd |Ho et al] (|2 004), 
and its pol arimetcr is described by Marrone' (200(5) and 
Marrone & Rao (2008j ) . The most significant detection of 
CP at 230 GHz was made from the observations of March 
31, 2007. Archival data from two other observations — 
one at 230 GHz and another at 345 GHz — were ana- 
lyzed to confirm and extend the initial result. A typi- 



^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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TABLE 3 

Double-Sideband Flux and (Fractional) Polarization for Sgr A* 



Date 


Freq. [GHz] 


[Jy] 


Q/I 


U/I 


V/I 


2005 June 6 


343.0 


3.17± 0.02 


0.027 ± 0.003 


-0.049 ± 0.003 


-0.016 ±0.003 


2006 July 17 


226.9 


3.88 ±0.02 


0.020 ± 0.003 


-0.053 ± 0.003 


-0.011 ± 0.003 


2007 March 31 


226.9 


3.52 ±0.01 


-0.061 ± 0.001 


-0.032 ±0.001 


-0.012 ±0.001 


2008 May 30'" 


226.9 


3.87 ± 0.02 


0.044 ± 0.002 


-0.032 ±0.002 


-0.011 ± 0.004 



Statistical error only. Absolute calibration precision is typically 10-20%. 



Sgr A* was observed for one hour prior to the 1924-292 test observation. 



cal Sgr A* polarimetry track consisted of ^ 12 hours of 
observations, of which ~ 6 corresponded to continuous 
monitoring of Sgr A*. The quasar 1733-130 was used as 
the gain cahbrator of Sgr A*, while quasars 3C273 and 
3C279 were used to cahbrate for instrumental polariza- 
tion (leakage). In addition, quasars such as 1337-129 and 
3C286 were commonly included in the tracks, as well as 
a solar system object (e.g.. Titan) for flux density cali- 
bration. 



M 




RA offset (arcsec) 

Fig. 2.— Total and CP emission of Sgr A* (March 31, 2007): 
red contours correspond to Stokes / flux density while blue con- 
tours correspond to Stokes V flux density. Solid and dashed 
contours indicate positive and negative flux density, respectively. 
Stokes / contour spacings are —6,6,12,25,50,100 X aj, where 
(Tj = 19.8 mjy/beam is the rms noise in the / map. The con- 
tour levels for the Stokes V map are -11,-7,-4,-2,2,4 X av, 
where ay = 2.9 mJy/beam is the rms noise in the V map. 



To test the accuracy of our CP measurement, we made 
a special observation of the quasar 1924-292 at 230 GHz 
on May 30, 2008, along with another short-track mea- 
surement of Sgr A*. 1924-292 has nearly the same dec- 
lination as Sgr A*, so both have the same AZ-EL track 
and, hence, position angle dependence with time. This 
track also contained a one-hour segment on Sgr A*. The 
gain calibrator for 1924-292 was 1911-201, and the polar- 
ization calibrators were 3C273 and 3C279. Every effort 
was made to keep the observational circumstances iden- 
tical to those of previous tracks (e.g., same polarization 
sampling cycle). 

The SMA polarimeter allows precise gain calibration 
for the RR and LL visibilities. A single quarter-wave 



plate is located in the beam of each telescope and ro- 
tates to produce either right CP or LCP. The polariza- 
tion of each antenna is modulated according to orthogo- 
nal Walsh functions of period 16 to efficiently sample all 
four polarization combinations on each baseline. Quasi- 
simultaneous polarization data are generated in post- 
processing by averaging over the switching cycle. The 
use of the same waveplate and feed for both R and L 
polarization states eliminates some of the uncertainties 
inherent in the use of dual-feed receivers, particularly dif- 
ferential phase variations between the polarization states. 
Measurement of LP (Stokes Q and U) relies on precise 
determination of the leakage of each polarization state by 
the crosshanded polarization, which was performed using 
linearly polarized bright point sources (quasars 3C279 or 
3C273) observed over a large range of parallactic angle. 
To first order, these lea kages do not affect the measure 



ment of Stokes V (e.g. , IMarrone 2006 



2001 



Thompson et al 
ative calibration 



and Section |4.1[ ) , leaving the re' 
of the RR and LL visibilities as the primary calibration 
challenge for measurement of Stokes V . 

3. RESULTS 

The polarization measurements (Stokes /, Q, U , and 
V) of Sgr A* for four different epochs included in our 
analysis are shown in Table [3j The flux densities for 
each Stokes parameter /, Q, U , and V were obtained 
by fitting a point-source model to the visibilities, fixing 
the position to the phase center of the observations, the 
location of Sgr A*. The polarization Stokes parameters 
are listed as a fraction of the total flux density. For four 
different epochs — including different frequencies, differ- 
ent Mw-space sampling, and a time span of three years — 
Sgr A* is shown consistently to be circularly polarized at 
the > 1% level. The sense of CP (negative V) persists 
throughout. 

Figure [2] shows the contour maps for both Stokes / 
and V for the whole track of March 31, 2007. In this 
case, the measured CP flux density is —41 ± 3 mJy (in 
the standard lAU sign convention, where a negative sign 
indicates LCP), corresponding to ~ —1.2% of the total 
flux density. The background rms noise of the image is 
3.1 mJy beam~^, consistent with the statistical error of 
±3 mJy obtained from the visibility fit. Figure |3] shows 
the compilation of published CP observations from the 
lower panel of Figure [l] along with our new measure- 
ments. Although there are many unobserved frequencies, 
the polarization appears to increase monotonically with 
frequency while retaining the same handedness through- 
out. The CP spectrum scales approximately as j/0-35±o.03 
across the range of frequencies with detections. Subdi- 
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Fig. 3. — Fractional CP in Sgr A* (LCP in all cases) from radio 
to submillimeter frequencies (lower panel in Figure [Tj , including 
the new SMA data (see Table Isll. The percentage ofTjP increases 
with frequency as ^ j/O.SSiO.OT"^ 

viding^he track into four segments, (see Figure [4] and 
Table 41, we see no statistically significant variability in 
the CP^ with the largest change corresponding to a 2-a 
difference. We limit fractional variation in the CP, i.e., 
AV/V, at the 2-a level of significance to 40%, which is 
comparable t o the fraction al changes in CP at 1—15 GHz 
measured by [Bower et al . | pO 02). There is a tantalizing 
correlation at the 2-<t level of significance between the 
circularly and linearly polarized flux densities. However, 
with only four data points, so significant conclusion can 
be drawn. 

Figure [2] shows a slight offset between the peak posi- 
tions of V and /. This offset, 0.18 arcseconds, is consis- 
tent with what is expected for an SNR of ^14. The offset 
is observed in all of the Sgr A* polarimetric tracks at all 
frequencies, and the orientation and angular amplitude 
of the offset vary with time within a single track. This 
shift between the peaks of the Stokes V and Stokes / 
images suggests imperfection in our CP calibration pro- 
cedure. The offset is barely noticeable for maps derived 
from complete tracks (e.g.. Figure [2]), but it becomes sig- 
nificantly larger when the track is split into consecutive 
time intervals (see Figure Inland Table [4]). The fact that 
the offset grows when reducing the SNR suggests that 
the wandering of the Stokes V map around the Stokes / 
may be in part due to a thermal noise effect. However, 
extensive tests (Section [4| suggest that the offset is not 
dominated by noise. 

Table [5] shows the polarization flux densities for all 
the test quasars observed on March 31, 2007. All point- 
source flux densities are obtained from the visibility fit- 
ting process. Quasars 3C273, 3C279, 3C286, 1337-129, 
and 1733-130 show CP fractions consistent with zero in 
both sidebands. The deviations from zero are of the same 
sign for all objects, which suggests a systematic error that 
would fractionally apply to Sgr A* as well. However, the 
0.1% weighted-average magnitude (right CP) of this sys- 
tematic error is negligibly small, an order of magnitude 
smaller than the observed CP in Sgr A*. 

The circular polarization of quasar 1924-292 (last line 
of Table [S]) was also measured in a separate track on 
May 30, 2D08. The V flux determined by a point-source 
fit to the visibilities, constrained to lie at the pointing 



center, is consistent with zero, as for the other sources. 
However, within a region comparable in size to the syn- 
thesized beam, the Stokes V signal varies between +5cr 
and —4(7, as shown in Figuref5] This antisymmetric pat- 
tern and other calibration effects are discussed in detail 
in Section mH 



4. ERROR INSPECTION AND POLARIZATION 
TESTS 

The consistency of the measurements in Table [3] over a 
period of three years and at multiple frequencies suggests 
a robust detection of CP. The most sensitive of the de- 
tections has a significance of greater than lOcr. Because 
our work presents the first measurements of CP with the 
SMA, we feel it is important to describe our calibration 
methods and system checks in detail. 

4.1. Polarization Calibration 

In an ideal interferometer with circularly polarized 
feeds, CP in a target source appears as a difference be- 
tween the parallel-hand visibilities. Complications may 
arise due to imperfection in the interferometer response 
and calibration. The four polarized visibilitie s in a circu- 
larly polari zed interferometer can be written ( Thompson 
|2001j ): 



et al. 



VbR = 9Rag*Rb [0^1 + 

+ dRa {Vq - iVu) e^' 



Vv)+dRad]^, (Vi-Vv) 
'' + d*Rb (VQ+iV[/)e-2^*] 



Vrl = 9Rag*Li, [(Vq + tVu) + dRa {Vl - Vv) 

-dlb {Vi + Vv) - dRadl, (Vq - iVu) e^^*] 

Vlr - gLa9*Rb [(^Q - ^^u) e^'^ - dia (V/ + Vy) 
+d*R, (Vi - Vv) - dLad*R, {Vq + iVu) 

VlL = 9Lagh [(Vl - Vy) + dLa^L (V/ + Vy) 

~dLa (Vq + iVu) e-2# _ (Vq - iVu) e^**] , 

where gRa and gLa (dRa and dia) are the right and left 
circular feed gains (polarization leakage terms) of an- 
tenna a, and </> is the parallactic angle of the feed. Solving 
the linear system of equations above for the Stokes / and 
the Stokes V visibilities, respectively, gives 



Vr 



Vrr 



(1 + dLadl,) 



'2gRag*Ri, (1 + dRadLa)(l + d*j^Jl^) 
VlL (1 + dRad*R^) 



2gLaglb(i^ 

Vrl 


- dRadLa)(l + d*j^^dl^) 
(dLa ~ dflfc) 


'^gRaglb (1 - 

Vlr 


^ dRadLa)(l + d*j^f^dl^) 
(d*Lb - dRa) 



(1) 



and 



Vy 



'^9La9*Rb (1 + dRadLa)(^ + d*R^d*^^) 
(1 - dLadlb) 



V 



RR 



V 



'^gRag*Rb (1 + dRadLa)(l + d*Rbdl^) 



LL 



(1 



dRad*Rb) 



Vrl 


- dRadLa)(l + d*j^Jll,) 
(dLa + dflj 


"^gRaglb (1 - 

Vlr 


1- dRadLa)(l + d*R^dl^) 
(dRa + dl^) 



'^9Rb9La (1 + dRadLa)(l + d*j^Jl^) 



(2) 
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TABLE 4 

Circular Polarization in Subsections of March 31, 2008 Track 



UT hour range 




V/I 




A6»v 
(arcsec) 


SNRv 


12.4 -14.2 


(- 


1.3 ±0.2) X 10" 


2 


1.52 ± 0.27 


6.5 


14.2 -15.7 


(- 


1.1 ± 0.1) X IQ- 


2 


0.78 ± 0.21 


7.8 


15.7 -17.2 


(- 


1.6 ±0.1) X IQ- 


2 


0.77 ± 0.10 


10.5 


17.2 -18.9 


(- 


1.6 ±0.2) X 10- 


2 


0.31 ± 0.18 


9.1 




10 5 -5 -10 10 5 -5 



RA offset (arcsec) RA offset (arcsec) 




10 5 -5 -10 10 5 -5 -10 



RA offset (arcsec) RA offset (arcsec) 



Fig. 4. — Contour images of Sgr A* made from four segments of the track on March 31, 2007. Red: Stokes / in Jy/beam. Blue: Stokes 
V in Jy/beam. The color coding and contour levels are the same as in Figure [2] Each of the four panels corresponds to one of the four 
time intervals the track was divided into: UT=13.4, 14.9, 16.4, and 17.9 hours. 



To first order in the leakages d, and ignoring terms pro- 
portional to Vlr d and Vrl d (small LP as well as small 
leakages), we have 



V, 



and 



=i ^ I VRR/{gRag*Rb) - "^LL/igLaglb)] , 



(3) 



(4) 



and thus Stokes / and V are independent of the leakages 
to first order. 



The MIRIAD reduction package ( jSault et al.||1995 l 
uses these first-order equations when solving tor the po- 
larized leakages, ignoring second-order terms in the leak- 
ages d and LP fraction. These terms contribute a sys- 
tematic error in Stokes V of the form Id^ and md, for 
LP fraction m. For the track of March 31, 2007, all leak- 
age terms are of order ~ 10~^ (both real and imaginary 



parts) , introducing a fractional contribution from Stokes 
/ of ^ 10^^. Likewise, the terms of order md are of 
order 10~^ and are unable to explain the observed CP 
fraction of 10~^. Another uncertainty arises from the 
complex terms dependent on the parallactic angle in the 
form of a phase term in the full gain equations. A system- 
atic phase effect could, in principle, explain offsets in the 
image plane; however, the low flux density contribution 
from these offset components (~ md) makes this possi- 
bility difficult to reconcile with the evident displacement 
of the entire point-source flux density in Figure |4] 

4.2. Gain Calibration 

Errors in the relative calibration of the L and R gains 
of the antennas affect the precise cancellation of terms 
in eq. [4] and are the most natural way to generate false 
CP. The time variation of the antenna gains is moni- 
tored using 1733-130, which, at a flux density of ~1 Jy, 
is much weaker than Sgr A*. In order to maximize the 



CIRCULAR POLARIZATION OF SCR A* 



7 



TABLE 5 

Circular Polarization'' for Test Quasars on March 31, 2007 



Source 




LSB 






USB 




I [Jy] 


V/I 




I [Jy] 


V/I 




30273*^ 


15.40 ± 0.03 


(1.2 ± 1.2) X 10" 


'S 


15.05 ±0.02 


(1.1 ± 1.3) X 10- 


3 


3C279 


12.92 ±0.02 


(1.5 ± 1.4) X 10- 


'3 


12.88 ±0.02 


(1.3 ± 1.4) X 10- 


3 


3C286 


0.49 ±0.01 


(7.9 ± 15.2) X 10- 


-3 


0.46 ±0.01 


(14.8 ± 17.9) X 10- 


3 


1337-129 


6.92 ±0.03 


(2.1 ± 3.0) X 10" 


-3 


6.89 ±0.03 


(2.5 ± 3.4) X 10- 


3 


1733-130'^ 


1.48 ±0.04 


(0.8 ±0.4) X 10- 


-3 


1.46 ±0.04 


(1.0 ±0.5) X 10- 


3 


1924-292^ 


7.00 ± 0.01 


(-0.2 ± 1.1) X ID- 


-3 


6.95 ±0.01 


(-0.1 ± 1.2) X 10- 


3 



^ Stokes V fluxes were calculated by fitting the visibility data to visibilities corresponding to a 
point source located at the interferometric (phase) center of the map. 

^ 3C273 was calibrated in polarization and gain with 3C279. In contrast, quasars 3C279, 3C286, 
and 1337-129 were calibrated using 3C273 as a gain and a polarization calibrator. 
^ Quasar 1773-130 is the gain calibrator for Sgr A*. The average gain curve was derived using 
both LL and RR visibilities, then applied to the 1733-130 data before measuring the CP. 
Shown here for comparison, quasar 1924-292 was observed on the night of May 30, 2008. 




RA offset (arcsec) 

Fig. 5.— Total and CP emission of 1924-292. Color cod- 
ing is identical to Figure [2] Stokes / contour spacings are 
12,25,50,100,200 x aj, where aj = 14.7 mjy/beam is the rms 
noise in the / map. The contour levels for the Stokes V map are 
—5, —3, —2, 2, 3 X ay where ay = 5.85 mjy/bcam is the rms noise 
in the V map. Within the angular extent of the Stokes / point 
source, two distinct Stokes V are identifiable. The negative peak 
corresponds to approximately —30 mJy/beam (i.e., —0.4% of the 
Stokes I peak), and the positive one corresponds to approximately 
+18 mJy/beam (i.e., 0.25% of the Stokes / peak). 

signal-to-noise of the gain measurements, we choose to 
average the amplitude of the LL and RR visibilities to- 
gether and determine a single gain curve per antenna 
rather than separately measuring and qr. If there is 
an imbalance in the response of the two "feeds," which 
are really the same receiver looking through the same 
quarter-wave plate, this will introduce false CP. 

The first test of our calibration method is to examine 
the calibrated 1733-130 data for CP. In Table[5]we report 
the polarization measured in 1733-130 after applying the 
same gains that are used with the Sgr A* data. Al- 
though we make no attempt to correct the difference in 
the L and R gains in our calibration procedure, we end 
up with little polarization in this source (<0.1%), with 
a fractional V sensitivity of ~ 5 x 10~^. This suggests 



a systematic polarization limit of 10^'^ for this calibra- 
tion method. Performing the calibration with separate 
gains for L and R yields nearly identical results for the 
polarization fraction of Sgr A* but with slightly higher 
noise. 

A second test of this calibration method is to examine 
the V/I measured for other quasars in the track. We 
expect no significant CP in quasars at these wavelengths. 
Figure [6] and Table |5] show that no CP is found in the 
quasars observed in the March 31, 2007, track. Fractional 
polarization limits are comparable to those obtained on 
1733-130, except for 3C286, which is significantly fainter 
than the others and thus has a higher level of noise in 
V/L 

4.3. Circular Polarization of 1924-292 

To guard against poorly understood instrumental ef- 
fects that preferentially affect, for example, sources at 
low elevation or those with poor parallactic angle cover- 
age, we undertook a supplementary test using the quasar 
1924-292. The declination of this source is within half a 
degree of that of Sgr A*, ensuring that it follows the 
same path on the sky. The observation scheme was de- 
signed to precisely mimic that of the primary track on 
Sgr A*, with similar hour angle coverage and temporal 
sampling, identical polarization modulation, and a gain 
calibrator (1922-201) at a similar distance from the tar- 
get (9° north compared to 16° for 1733-130 and Sgr A*). 
Because 1924-292 was brighter at the time of these ob- 
servations, 7 Jy, than Sgr A* was during the other ob- 
servations (3-4 Jy), and because of the long integration 
time in the track, this test also provides a more sensi- 
tive limit on instrumental CP than the that provided by 
observations of quasars on March 31, 2007. 

The results of this test are shown in Figure [5] and the 
last line of Table [5l Fitting the visibilities to a point 
source, as was doneior other quasars and Sgr A*, we find 
V/I = 1 - 2 X 10"'', with an uncertainty of 10"^. This 
is a factor of more than 10 lower than the CP fraction 
observed in Sgr A*, including that measured in the same 
track in a short observation (last line. Table [Sl. That 
the polarization of 1924-292 was not found at the same 
time that the previously observed CP level in Sgr A* was 
reconfirmed enhances our confidence in our result. 

Figure [5] does show significant CP in an antisymmet- 
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Fig. 6.— Quasars 3C273, 3C279, 3C286, and 1337-129 observed during the track of March 31, 2007, after gain calibration using 3C273 
as cahbrator. Each image shows contour plots in Stokes / and V in the same way as Figure [2] For each panel, the contours in Stokes V 
are -2 and 2 times the rms noise of the image, and all four of them have maxima and minima between —"iuy and +3(tv. The intensity 
contours in Stokes / are (-2,2,5,10,20,40,80,160,320,640) X ai for 3C273 and 3C279. For 3C286 and 1337-129, the contour levels are 
(-2,2,5,10,20,40) and (-2,2,5,10,20,40,80, 160) X ct/, respectively 



ric pattern around the pointing center, with a peak V jl 
values of -1-2.5 and —4.2 x 10"'^. The negative peak is 
approximately five times the image rms. The antisym- 
metric pattern is clearly distinguishable from the offset 
CP detections in Sgr A*, which do not show peaks of 
opposite sign one resolution element away from the main 
peak. We conclude that this is a signature of an unre- 
lated calibration problem. This calibration could mask 
a peak in the V map that is offset from the pointing 
center as seen in Sgr A* in Figure 111 with the allowed 
amplitude being approximately the difference in the pos- 
itive and negative peaks, or —1.7 x 10""^. We take this 
as a limit on our uncorrected instrumental polarization, 
y// < 2 X 10-3. 

The LP of 1924-292 was measured to be 10.8% dur- 
ing this track, larger than was observed in the Sgr A* 
tracks in this work. The increased LP makes this source 
slightly more susceptible to errors from ignored second- 
order terms from feed imperfections (Section |4.1[ ). Nev- 
ertheless, no significant CP is detected. 

4.4. Other Possible Sources of Error 

In the tests discussed thus far, we failed to find CP in 
any source other than Sgr A*, placing tight limits on the 
possible instrumental contribution to the observed signal. 
The previous section describes the most stringent test, 
which mimics everything about the Sgr A* observation 
and nevertheless finds insignificant CP in a bright test 



source. 

A possible distinction between the test sources and 
Sgr A* is the presence of extended emission in our Galac- 
tic Center. The test sources are nearly pointlike (with 
some low-level emission from jet knots in the case of 
3C279), while Sgr A* is embedded in dust and free- free 
emission. To avoid contamination from larger-scale emis- 
sion, we exclude baselines shorter than 20 kA from our 
visibility fits and find that the same exclusion in other 
sources does not change our results. We are not aware 
of a mechanism by which the extended emission could 
introduce false CP, but we have also attempted to ver- 
ify that no such effect is present. We examined a data 
set from July 17, 2006, which measures uv spacings of 
27-390 kA instead of the 5-50 kA that is typical of the 
other tracks. This track should be more immune to the 
larger scale emission, which is highly suppressed by the 
shorter fringe spacing. We see no significant difference 
in Table 13] between this track and the others. 

5. ASTROPHYSICAL SOURCES OF CIRCULAR 
POLARIZATION 

The polarization state of radiation emerging from an 
astrophysical source is governed by radiative transfer 
equations that, in their most complete form, incorporate 
emission and absorption in each Stokes parameter. This 
polarized radiative transfer equation must also include 
mixing between Stokes parameters that takes place in 
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birefringent media, such as a synchrotron plasma. The 
magnetic field that is fundamental to the synchrotron 
emission and self-absorption is also a source of birefrin- 
gence. Variations in the magnetic field along the line of 
sight, and the change in the coupling between Stokes pa- 
rameters with plasma temperature, complicate the pic- 
ture further. A fully consistent treatment requires ad- 
dressing all these effects simultaneously. 

5.1. Polarized Transfer in a Homogeneous Magnetic 

Field 

Sg r A*, a stratified , self-absorbed synchrotron source 



(see de Bruyn 1976 for the basics of these models) 
with frequency- and time- variable polarization, clearly 
requires a sophisticated radiative transfer model. The 
Farad ay rotation inferred from submillimeter observa- 
tions ( Marrone et al. 2007 1 has been assumed to be sep- 
arate from the submillimeter emission regions, although 
it may instead occur deep inside the source ( sometimes 
referred to as internal Faraday rotation, e.g., W ardle "fc| 
Homan 2003). The production of CP may similarly be 



a signature of conversion processes within the emission 
region . The complete PRT equation is (e.g., Sazonov 
1969| 







( 111 






VQ 


dl 




i 






\Vv 




-Hj —Kq — Ky 
-kq -Hi -Ky 

+Hy —Kj —K*Q 

'Kv +Kq —Hi 

where the cartesian axes in the plane of the sky (ei and 
62) — which determine the precise definition of Stokes Q 
and U — are oriented such that ei is aligned with the local 
direction of the magnetic field. The absence of transport 
coefficients related to Stokes U in eq. [5] is a direct conse- 
quence of this convenient choice of coordinates, and thus 
it is a valid expression only locally. 

Eq.[5]is the vectorial generalization of the standard ra- 
diative transfer equation dl^/dl = ju — a^I^. The rjA 
(with A = I,Q, U, V) are the emission coefficients for 
each Stokes parameter. The absorption coefficients ka 
(one for each Stokes parameter) are related to the tja 
through detailed balance. The antisymmetric portion of 
the matrix comprises the Faraday transport coefficients 
K*A- They couple Q, U, and V, permitting exchanges 
of LP and CP along the path of radiation. This ro- 
tation in three-dimensional polarization space (Q, U, V) 
has b een referred to as the generalized Faraday rotation 
effect (Melrose 1997). The general Faraday effect reduces 
to pure l^'araday rotation (i.e., rotation in the Q-U plane) 
in the nonrelativistic limit and to pure Faraday conver- 
sion (i.e., rotation in a plane perpendicular to the Q-U 
plane) in the ultrarelativistic limit but, in general, is a 
combination of both. 

The intrinsic emission term in eq. p](7;y) demonstrates 
the first potential source of CP. The intrinsic emission 
in Stokes V is weak: rjy is of order I/7 smaller than the 
other components and is generally ignored in synchrotron 
emitting plasma. On the other hand, predominantly CP 
emission can be expected in a cold c yclotron emitting 
plasma (e.g.. Landau fc Lifshitz||1975 1, in which the rel- 
ativistic approxima tions intrinsic to synchrotron spectra 
(e.g., beaming, see Rybicki & Lightjnan 1979) are not 



valid (see also |Mahadcvan et al.^^1996 ). ihe submillime- 



ter emission from Sgr A* is expected to arise at very small 
radius where the plasma is hot, and therefore emission 
from cold electrons is not expected to be relevant. 

The birefringence of the plasma may also transform LP 
to CP, as noted above, but this "Faraday conversion" 
is weak in a cold plasma. In this limit, the response 
of the medium depends exclusively on the plasma fre- 
quency, Vp = (nee^)^/^/(7rme)^/^, and the cyclotron fre- 
quency, = 6-6/ (27rmec). In the high frequency regime 
(i.e., vb/i^ <C 1, applicable to the submillimeter emission 
considered here given plausible magnetic field strengths), 
the rotation and conversion coefl ncients are (Melrose & 
McPhedran|[2005l |SwansMl|[T989l ) 

(6) 
(7) 



■■{2TT/c)iuliyB/'y'^)cosl 



sm 



where 6 is the angle between the magnetic field vector 
and the line of sight. The ratio of these equations shows 
Faraday conversion (Ky) to be a factor of vb/i^ (~ 10~'^ 
for our data) smaller than rotation, inconsistent with our 
observations. Equations for higher temperature plasmas 
are presented in section [5. 6| 



5.2. Magnetic Field Changes 



iHoman et al. (2009) list a series of scenarios in which 
CP can be produced in quasars. By incorporating line-of- 
sight changes in the magnetic field geometry, they offer 
possibilities not present in the homogenous-field radia- 
tive transfer equation (eq. [5]) . In a uniform field. Stokes 
Q needs to be converted into Stokes U before Stokes U 
can be converted into Stokes V. Changes in the field 
orientation along the line of sight break a fundamental 
assumption of eq.p] the azimuthal angle (</>) between the 
field direction in tne plane of the sky and the coordinate 
axis ei was defined to be zero. Changes in field direction 
will introduce a coupling of Q and V across the domains 
of magnetic orientation. The magnitude of the transfer 
coefficients also depends on the angle 9, introducing a 
small additional change in couplin g. These field orienta- 
tion chang es, whether stochastic ( Ruszk owski &: Begel 
man [| 2002 ) or ordered (as in the helical fields ot |Beckert 
20031) can also create significant CP. Both of these mech- 



amsms are capable of producing the CP sign coherence 
observed in Sgr A*, although a small net bias to the field 
direction is important for the stochastic field. 

5.3. Self-Absorbed Faraday Conversion 

The wavelength-dependent size of Sgr A* indicates 
the importance of synchrotron self- absorption in the ob- 
served SED. The source size and brightness are deter- 
mined by the variation of density, magnetic field, and 
electron energy distribution with radius/distance, al- 
though there are unbreakable degeneracies between these 
quantities even assuming power laws in radius with- 
out additional information or assumptions. A specific 
choice of jet structure can determine the structure in 
density / field / energy and match both size and flux (e.g 



Markoff et al. 



eke et al. 2005 



an 



20071 and even variability timing (|FaI 
Himilarly, the assumption of, e.g, 
equipartitlon magnetic field can specify the structure for 
an accretion flow while matching observations. 

Polarimetric observations add new information regard- 
ing the source structure, potentially reducing the need 
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for unverifiable assumptions. The structure of Sgr A* is 
likely to be more complex than is captured by power- 
law models; this is particularly true at submillimeter 
wavelengths where the emission originates near the black 
hole. Complex magneto-hydrodynamic simulatio ns and 
general relativistic radiative tra nsfer (e.g., , Huang et al. 
2008; Moscibrodzka ct al. 2009; Shcherbakov et al.||ij010| ' 
S!icherbakov fc Huang„2011^ will more faithfully repre- 
sent the true source properties near 345 GHz than few- 
parameter analytic models. Nevertheless, the smooth- 
ness of the average CP spectrum, the exclusive left- 
handedness of the polarization from 1.4 to 345 GHz, 
and the slow monotonic variation of the fractional po- 
larization suggest that much insight can be gained from 
a simple conceptual model incorpor ating the ideas of self- 
absorbed synchrotron sources (e.g., |de Br{I^|1976[ ). 

In the case of self-absorbed synchrotron sources, the 
specific intensity near frequency v is dominated by the 
emission of electrons within a narrow energy range 
around 7rad- An electron with Lorentz factor 7i.ad 
emits most of its radiative power at a frequency Vc ~ 
^"^^^eB sinO / {2'nmf.c) = ^"^^^vba.- The brightness tem- 
perature associated with the emitting elec trons is re 
lated to this effective Lorentz factor by 



1970 



Rybicki & Lightman 1979) ksTb = a{p)jra.d'm.eC 



le.g., 



Pacholczyk 
2" 



(see, for example, Loeb & 



where a[p) is a coetticient ot order unity. The brightness 
temperature can be obtained directly from observations: 
Tb = 1.22 X lO^^S^v-^ej^iv) K; where is the ob- 
served flux density in Jy, is the source angular radius 
in milliarcseconds, and ly is in GHz. Assuming power-law 
dependencies of the form oc i/^™ and 6's(j^)_oc 
we obtain 7i.ad oc i,-"i-2+2n 
[Waxman 2007 ). 

5.4. Self-Absorbed, Stratified Sources 

In what follows, our main approximation is to assume 
that at each frequency, the synchrotron-emitting elec- 
trons, with energies 7radj are also responsible for Fara- 
day conversion. We also assume that Faraday rotation 
and conversion can be considered to act in an alternating 
manner rather than attempting to treat the full, compli- 
cated PRT equation. The radiative transfer near the 
r = 1 surface is particularly complicated because all ma- 
trix elements can be important. By assuming sequential 
effects, we can operator-split the differential equation. In 
this case, and for short propagation distances, the pro- 
duction of Stokes V is simply an angle of rotation in the 
Poincaire sphere corresponding to the cumulative phase 
shift betw een the polarization mode s of the plasma near 
radius rg ( Kennett fc Melro"se||1998 ). This can be calcu- 
lated as an integral ot Kq, 



A')p{ro,iy) 



sin 9 9fc dr 



(8) 



where we have corrected the cold-plasma Faraday conver- 
sion rate (eq. [t]) by a factor 9fc (see Section 5.6 below) 
that takes into account the effect of a power-law distri- 
bution of relativistic electrons on Faraday conversion. 

Followi ng the treatment of Farad ay conversion pre- 
sented by Kennett & Melrose (1998), we take the finite- 
temperature correction to be proportional to the local 
mean Lorentz factor of the electron distribution, and we 



same approximation 9fr, y log (7) 7(7); e.g., Sazonov 
1969| iQuataert fc Gruzinov||2000i , 

eFc(^^,Te,0,i3)^ (7(r)) ^7rad • (9) 

Assuming conversion in a narrow range in r, we write 
AV'(ro,i^)oc / ne{r) B (r)"^ ain^ jrad{r)iy^^dr 

Jro (10) 

ne{ro)B{ro)'^ jra,d{ro)v^^ro . 

We can obtain Aijj as a function of frequency by using 
the fact that for each observed frequency vq, there is a 
corresponding radius tq. If we assume, as well, that the 
electron density and magnetic field have radial profiles 
of power-law form, rig ~ r~^ and B ^ r~", the Faraday 
conversion phase shift is 



-0-2Q+1 



i+2n-5 



(11) 

Power-law relations of the form S^ oc v — for the 
flux density — and 6s{v) oc j/^" — for the source angu- 
lar size — are a na tural conseque nce of stratified syn- 
chrotron sources (de Bruyn 1976). In this simple 
model, the size-frequency relation is obtained as follows. 
First, we take the synchrotron op tical depth associated 
with a narrow shell of radius r ( [Rybicki fc Lightman 
19791 : T^(r) = a^Ar - ne{r)B'^P+'^1/'\r)v-^P+'^1/'^r 

a{p+2)/2 



-(p+4)/2 



where p is the electron power- 



law index. Hence, for a given frequency Vq, the r = 1- 



surface occurs at r. 





-(p+2)/(2/3-|-a(p-|-2)-2) 



-a(p-(-2)/2+l -(p-h4)/2 



1 or To 



^(13-50-3a-2ap+2p)/(; 

tral indices are 



Similarly, th e flux density can 



be shown to be (see de Bruyn 1976 for details) Sv 



The resulting spec- 



ie 



2/3- 
13 



a{p 

-5/3- 



h2) 
3a 



2 

2ap - 



and 



2p 



(12) 



2 - 2/3 - 2a - 



Therefore, a and j3 can be solved for as functions of m 
and n: 



-3 — 2m + 5n — 2p — mp + 2np 



2m — An 



and 



(13) 



which gives (dropping the subscript in vq for clarity) 



All; 



with l^{2n-m- 2){p - 1) 
(p-l)(-a + /3-l) 



(14) 



a(p + 2) + 2(/3 - 1) 



The indic es a and /3 are very sensitive to n, which is 



close to 
1.3 ± 



1, oE 



Falcke et al. 



( |2009[ ) report a value of n = 
tained from observations between 22 GHz 
and 230 GHz, after correcting for interstellar scattering. 
These results are very sensitive to the scattering model. 
Another approach, used by Shcn (2006), is to minimize 
the dependence on the scattering law by relying on the 
shortest wavelength data. For measurements at 43 and 
86 GHz, their result is n = 1.09 ±0.3. An even shallower 



further approximate this with our value of 7rad (in the slope can be obtained using the results from Doeleman] 
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et al. (2008). From the VLBI visibilities, the angular 
size of the source, assuming a circular gaussian profile, 
is 37 fias, which is smaller than the last photon orbit 
for a Schwarzschild black hole (~ 52 ^ta s). If the source 
geometry is assumed to be an annulus (Doeleman et al. 
[2008) ), the mean diameter is 58 /ias. Using this angular 
oiameter at 230 GHz gives a value of rt ~ 0. 8. In par- 
ticular, for this value of n, with m — —0.43 (An et al. 



2005 



we o 



and assuming an electron power-law index of 2.4 
□tain a ~ 1.2 and f3 ~ 1.4. These values are close to 
those corresponding to spherical accretion and equiparti- 
tion of energy (/3 = 3/2 and a = (/3 + l)/2 = 5/4). Most 
importantly, these indices imply I ~ 0.1, i.e., the CP pro- 
duction rate is almost independent of frequency, which is 
consistent with observations (ranging from 0.5% to 1.5%) 
across a range in frequency of almost two orders of mag- 
nitude. The main point here is that a mild increase of 
CP with frequency is possible. Although this is contrary 
to the decreasing efficiency of Faraday conversion with 
frequency that is expected for a uniform medium, it is a 
natural consequence of stratified emission. 

5.5. Magnetic Field Structure and Polarization 
Stability 

The amount of CP generated from LP is proportional 
to this phase shift, i.e., rric ^ mi Ai/j for small Aip. The 
weak dependence of Aip on frequency implies that CP 
production in different layers of the source behaves sim- 
ilarly regardless of location respect to the SBMH. The 
value of Alp should remain small to avoid conversion 
of CP back to LP at high freque ncies. Starting with 
the ~ 10% (e.g., Pacholczyk||1970 ) intrinsic polarization 
fraction of the optically thick synchrotron, achieving CP 
fractions of 0.1% to 1% requires Aip ~ 0.01-0.1. Such 
small values are far from reversing the sign of CP, con- 
sistent with the fixed sign of CP at all frequencies, as 
long as the orientation of the orientation magnetic field 
does not fluctuate randomly. A variable orientation of 
the magnetic field with radius — i.e., field reversals — is 
inconsistent with the apparent self-similarity of the CP 
spectrum. If this self-similarity persists for more than 
two orders of magnitude in frequency (1 to 345 GHz), the 
orientation of the magnetic field must be highly coherent 
over spatial scales spanning nearly two orders of magni- 
tude (since n « 1), roughly between IRs and 100i?5. 

In plane-of-the-sky coordinates aligned with Bg^y (to- 
tal field equals plane-of-the-sky field Bg^y plus line-of- 
sight field Bios), eq.[5]shows that for intrinsic synchrotron 
emission only in / and Q, Faraday conversion proceeds 
as Q U —T' V, i.e., conversion is driven by Faraday ro - 
tation within a narrow synchrotron shell ( Beckert||2003 ) . 
The sign of U generated by this mechanism depends on 
the sign of the Faraday rotation coefficient K,y (eq. [6]), 
which can change depending on the angle the magnetic 
field makes with respect to the line of sight 9. On the 
other hand, the sign of the Faraday conversion coefficient 
Kq (eq.lTl) does not change with 9. This implies that the 
sign of Stokes V generated by conversion is exclusively 
determined by the sign of Stokes U generated from Q 
by local Faraday rotation, and is ultimately dependent 
on the local value of cos 9. Therefore, a constant sign in 
Stokes V needs a constant orientation of Bjos. This re- 
quirement is a lready hinted at by th e observed RM prop- 
erties ( jMarrone et al.|[2006a 2007). On the other hand. 



the orientation of Bgky is free to change in time without 
affecting the value of V. This is a consequence of the se- 
quence Q U V given by eq. |5] where the first step 
is performed by Faraday rotation depending only on Bjos 
and the fact that Stokes V does not change with coor- 
dinate transformations in the plane of the sky. Changes 
of the orientation of Bgky in time will only change the 
observer's definition of Stokes Q and U. Indeed, this is 
likely to be the exp lanation for the observ ed wandering 
values of Q and U (Marrone et al. 2006b I, while Stokes 
V appears to be stable. 

5.5.1. Relation to Observed Faraday Rotation Rates 

A key assumption in the self-absorbed approximation 
above is that synchrotron emission, Faraday conversion, 
and Faraday rotation (which drives conversion) all take 
place within t he same narrow shell of plasma. Previ- 
ous work (e.g., Marrone et al. 2007) has assumed that 
Faraday rotation primarily occurs tar from the emission 
region at a given frequency. This "Faraday screen" ap- 
proximation may still operate in regions of the accretion 
envelope where the material is too optically thin to con- 
tribute any emission or absorption to the incoming radi- 
ation while still being dense and magnet ized enough to 
produce Faraday rotation. For such cases, [Marrone et al.| 
( j2006 a) calculate the cumulative Faraday rotation acting 
upon polarized radiation emitted at radius ro (analogous 
to eq. [si as 



Ax(ro,J^) = 



2'KvivB 



cos 9 Ofr dr , 



(15) 



(except for the 8fr term, which is taken to be 1). The 
total emission (Stokes /, V, and total LP) is assumed 
to come from a region interior to rg and to be essen- 
tially unaffected as it emerges from the source, except 
for the Faraday rotation of the electric vector position 
angle (EVPA). 

To be consistent with our self-absorbed Faraday con- 
version approach, Faraday rotation should act both in- 
ternally (at r ^ 1, where it drives conversion) as well 
as externally (at r ^ 1, where the Faraday screen ap- 
proximation is valid). To satisfy this requirement, the 
Faraday rotation scale lengths should extend far beyond 
the thickness of the emission/conversion shell. To con- 
firm this, we define the differential Faraday rotati on a nd 
conversion depths (differential forms of eqs.jsjand 15 1 as 



FR 



dr 



oc r-(^+")eFR(r)i.-2 « ^-(ff+a) ln7rad(r) 

^ — J (r) 



7rad(r) 



and 



(16) 



cJtfc 
dr 



oc r-(^+2")eFc(r)^.-^« r-(^+2")7,,d(r)rg" 



(Section 5.3) is also 



(17) 



respectively, where 7i.ad r 

a power law of r, with S small. It can be read- 
ily seen that dr-pc/dr is significantly steeper in radius 
than dr-pn/dr. For comparison, the synchrotron opti- 
cal depth can be written in differe ntial form as dr/dr ^ 
^-,3-a(p+2)/2j,-(p+4)/2 (ggg Scction [sisi . For values oi p 
between 2 and 3, dr/dr can be as steep as or steeper than 
drpc/dr, and thus also steeper than dTpsi/dr, confirming 
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that significant Farad ay rotation takes pla ce outside of 
the T — 1 surface (see Jones fc Odell|[T977 ). This means 
that Faraday rotation is still at work when Faraday con- 
version has stopped being effective. In other words, the 
source has only a limited spatial range to produce Stokes 
V — driven by the sequence Q ^ U ^ V — before the 
conversion rate decreases to negligible values. Therefore, 
while Faraday conversion is roughly local to the shell of 
radius tq, Faraday rotation is not. 

A consequence of having significant Faraday rotation 
within and immediately outside the photosphere is a 
frequency-dependent rotation measure. The highest fre- 
quencies probe deeper layers in the stratified source 
and therefore integrate Faraday rotation contributions 
from the shells that define lower-frequency photospheres. 
There are, as of yet, few observational constraints on this 
possibility. The RM has been determined through the 
measurement of a change in y acr oss a narrow range 
in frequency (Marrone et al. 20071, and less securely 



through the comparison of position angles averaged over 
sever al observations at mo re widely separated frequen- 
cies ( Macquart et al.||2006 ). A demonstration of the fre- 
quency dependence of the KM could be obtained through 
position angle measurements at multiple pairs of closely 
spaced frequencies, which should be possible very soon 
with improved millimeter interferometers and ALMA. 

Although the Faraday screen approximation is reason- 
able for a stratified synchrotron source, it can be difficult 
to obtain an accurate RM observationally, due to the 
source's layered geometry. Assuming the cold, optically 
thin Faraday rotation relatio n x(j^) = Xo + c^/iy^KM 



(e.g., Marrone et al. 2006a eq. 1), RM is observa- 
tionally obtained by computing the quantity dx/d\^ = 
dx/ d{c/v'^) ~ Ax/A(c^/i^^), using two neighboring fre- 
quencies and their respective observed EVPAs. The im- 
plied assumption, in addition to the validity of the screen, 
is that the emission at the two observed frequencies is 
originated in the same region and with the same initial 
EVPA xo- However, in the layered scenario, if these two 
frequencies correspond to two shells that do not over- 
lap spatially, they will produce emission that will travel 
through two different Faraday screens. In that case, RM 
would be a function of the two frequencies observed and 
not a constant value. Thus far, observations suggest that 
this is not the case and th at RM is indeed a c onstant, 
even when xo varies in time ( Marrone et al. 2007 1 . Simul- 
taneous observations at more than two frequencies will 
help resolve the question of whether Faraday rotation in 
Sgr A* is proportional to v~'^ or not. 

One unaddressed problem in the discussion above is 
that there apparently is not enough observed LP at low 
frequencies, (e.g., 4 and 8 GHz, where CP dominates) 
to generate the observed Stokes V via the Faraday con- 
version process. However, synchrotron flux density is 
highly polarized for ordered magnetic fi elds, even in the 



optically thick limit (Pacholczykl 1970 ). Therefore, it is 
highly likely that the observed LP suffers severe beam 
smoothing at these frequencies, at a point beyond the 
conversion region, if the Faraday-rotation region shows 
fluctuations (e.g., stochastic or turbulent in nature) at 
spatial scales smaller than the source size. For the sake 
of argument, consider partially polarized radiation with 
constant polarization angle xo entering a piecewise uni- 



form Faraday screen. Radiation will exit the screen with 
rms fluctuations in EVPA equal to ^ o-rm{c/ i')'^ , 
where ctrm contains the fluctuations in Ue, B, and 9. 
If the fluctuations in the outcoming Ax are Gaussian, 

9 2 

the observed LP will be {mi) = mifi e~ '^^^ where m;_o is 
the original LP fraction entering the screen. This occurs 
because, while Q and U are additive quantities, total LP 
is not, resulting in observed values of mi that are severely 
sensitive to EVPA changes in the plane of the sky. 

5.6. The Faraday Conversion at Finite Temperatures 

In an accretion flow, many of the physical quantities 
relevant for synchrotron emission are expected to in- 
crease inwardly. This is true for the electron density 
and the temperature T. In particular, the temperature of 
the plasma near the event horizon can achieve mildly to 
highly relativistic temperatures. It is evident then that 
the cold plasma approximation should not apply in the 
interior regions where Faraday conversion could operate. 
A more general expression for the Faraday transport co- 
effcients (eqs. |6]and[7]) is 



27r^'>s 



cose'eFR(i^,re,6',5) 



''Vcold 



QM^,T,,e,B) 



2 ? 



s\n'eQyc{v,T,,e,B) 



^4cold©Fc(i^,T„0,i?) 



(18) 



(19) 



where the functions OfrXi^, Te, B) and Q^ci^: ^e, B) 
are flnite-temperature corrections to the standard cold- 
plasma transport coefficients. By definition, 6fr,fc ^ 1 
when T — > and v ^ vb- 

In a cold plasma, the electrons oscillate in response 
to electromagnetic waves but are effectively stationary 
otherwise. When the electron energy distribution is non- 
negligible, the interaction of electron motions and inci- 
dent electromagnetic radiation is given by the Boltzmann 
equation in the presence of the Lorentz force: the Vlasov 
equation. The Vlasov equation describes the dielectric 
tensor of the plasma and thus its Faraday transport coef- 
ficien ts, which depend on the electr on energy distribution 
(e.g., [Melrose fc McPhedran||2005[ ). 

One of t he earliest app roaches to this problem origi- 
nated with'Sazonov ( 1969 ). For a power-law distribution 
of electron energies, the ratio of the Far aday transport 
coefficients (see also Jones &: Odell|[l977) is 



pi 



e 



pi 

FC 



C)Pl 



cold "^FR 

2(P+1) 



cold 



^3. 

/ mm 



(p + 2)(p-2) ln7„ 



2 VB^ 

I mm 



(p-2) 



(20) 



which is an increasing function of 7in n_(7min > 1 always 
under the condition vbIv 1. Eq. 



20 shows that even 



with ('*q/'*i/)j.qij ^ 1; for tiot enough plasmas, the Fara- 
day conversion coefficient Kq can become a significant 
fraction of Ky or even exceed it in extreme cases. 
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For electrons in a relativistic thermal (Maxwell- 
Juttner ) distribu t ion, t he plasma dielectric tensor is dif- 
ferent. Melrose (19971 derived an ultrarelativistic ap- 
proximation tor the response tensor of a magnetized ther- 
mal plasma. Ballantyne et al. ([2007) combined this re- 
sult with that ot a classical cold plasma to obtain an ap- 
proximate continuous function of th e polarization proper 
modes as a function of temperature. Shcherbakov ( 2008 1 
extended these results, providing consistent analytic ex- 
pressions valid for all temperatures. In this case, the 
ratio of the Faraday transport coefficients is 













therm 




cold 



©therm 
FC 

therm 



-V 



Did 



) 



+6- 



Kr 



\kFr ) 



g{T,9, ub/v) 



(21) 



where Kn {x) is the modified Bessel function of the sec- 
ond kind. Eq. [2l] smoothly reduces to the cold plasma 
limit when T — > 0, since the term in the square brackets 
becomes exactly 1 in this limit. At temperatures of lO^K 
(fcT/mgC^ = 1/6), the term in the square brackets takes 
a value of - 2.5, and at 2 x 10"K (kT/niec'^ = 34), it 

reaches ~ 1.3 x 10^. 

The factors / and g, above, introduced by |Shcherbakov 
(2008), take a value of 1 for vb/i^ — >■ and provide 
the necessary corrections for higher-order terms in vbIv- 
They can lead to significant modifications to the ra- 
tio of the Faraday transport coefhcients in the mod- 
erately relativisitic regime. For example, for intense 
magnetic fields, the cyclotron frequency can be writ- 
ten vb = 0.35 GHz (B/20G), and thus the ratio vsjv 
takes values of ^ 10"^ ^ at_submillimeter wavelengths. 
Following Shcherbakov] ( |2008[ ), at vbIv = 10^^ and 
T — 2 X 10^^ K with 9 = 7r/4, the correction factors 
become g ~ 0.93 and / ~ —0.003. The multipliers of the 
cold Faraday transport coefficients in eq. 21 can strongly 
change the expected behavior even at moderate temper- 
atures. This clearly demonstrates the importance of in- 
cluding these factors in numerical simulations of Sgr A* 
in order to properly handle the PRT problem. 

6. CONCLUSIONS 

We have reported the first detection of circularly po- 
larized radiation toward Sagittarius A* at submillimeter 
wavelengths. The detected CP fractions are at levels of 
1.2 ± 0.3% and 1.6 ± 0.3% for 230 GHz and 345 GHz, 
respectively, and are confirmed by observations at differ- 
ent epochs. We have tested the significance of this detec- 



tion by analyzing the gain calibration, spectral variabil- 
ity and time variability, and by observing other sources. 
CP is not found in any other source, with the most sen- 
sitive limits confining uncalibrated instrumental CP to 
less than 0.2%. 

Our measurements of CP, combined with previously re- 
ported measurements at lower frequencies, indicate that 
the polarization fraction rises with frequency. The sign 
(handedness) of the CP signal is the same for all de- 
tections at all frequencies over the period 1981 to 2007, 
as w as found in low- frequency measurements alone (e.g., 
[Bower et al.||2002[ ). The average CP fraction as a func- 
tion ol frequency is remarkably well characterized by a 
power law with i/0.35±o.03_ However, there have been no 
detections of CP between 15 and 230 GHz, so the CP 
spectrum may not be monotonic. 

The general trend of slowly increasing CP with fre- 
quency is consistent with self-absorbed Faraday conver- 
sion. We discuss a simple set of power-law scalings of 
B and Ue that allow Faraday conversion to operate with 
only weak frequency dependence, imposing no CP sign 
reversals as long as the magnetic field itself does not 
reverse over the range of radii providing centimeter to 
submillimeter emission. In contrast, more sophist icated 



models in the literature ( Huang et al.|2008 2009a) show 
frequent reversals in submillimeter CP and therefore are 
not consistent with the results presented here. 

If the CP originates from self-absorbed Faraday con- 
version, the corresponding Faraday rotation depth could 
vary with frequency, which would be contrary to the as- 
sumption of a cold Fa raday screen ( Macquart et al.|2d06 



Marrone et al. 2007). Polarization observations in the 
near future will test the frequency variation of the RM, 
and there already exists a clai m of an RM change be- 
tween 22 and 230 GHz (Yusef- Zadeh et all [2007 ). In 
general, the PRT in a self-absorbed synchrotron source 
should couple the variations in LP and CP, so sensi- 
tive observations of variability in the polarization of this 
source will be crucial to understanding the structure of 
the emission region and the source of CP. As with varia - 
tions in the RM (Shar ma et al.||2007| [Pang et al.||2011[ ), 
variations in CP will occur on a timescale set by the 
structure of the accretion flow and can therefore differ- 
entiate between models. 
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